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ABSTRACT 

The crustal deformation patterns associated with the 
earthquake cycle can depend strongly on the rheological 
properties of subcrustai material. Substantial deviations from 
the staple patterns for a uniformly elastic earth are expected 
when viscoelastic flow of subcrustai aaterial is considered. The 
detailed description pf the deformation pattern and in particular 
the surface displacements displacement rates, strains, aud 
strain rates depend on the structure and geometry of the aaterial 
near the seisaoge^ic zone. In the past fet[ years various 
viscoelastic models of crustal deformation have been publi id. 
These models differ in their predictions concerning the temporal 
and spatial patterns of crustal deformation. In some cases these 
differences are due to varying choices for the physical mechanism 
under study. In other cases^ however, the differences are the 
result of the details of the mathematical treatment or the choice 
of model parameters. We seek to resolve the origin of some of 
these differences by analyzing several different linear 
viscoelastic models with a common finite element computational 
teeijnigue. The models involve strike-slip faulting and include a 
thin channel asthenosphere model, a model with a varying 
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thickness lithosphere, and a model with a viscoelastic inclusion 
below the brittle slip plane. lh§ calculations reveal that the 
surface defoliation pattern is aost sensitive to the rheology of 
the material that lies below the slip plane in a volume whose 
extent is a few times the .fault depth. If this material is 
viscoelastic, the surface d,e formation pattern, resembles that of 
an elastic layer lying over a viscoelastic half-space. Uhen the 
thickness or breadth of the visaoelastic material is less than a 
few times the fault depth, then the surface deformation pattern 
is altered and g.eodetic measurements are potentially useful for 
studying the details of subsurface geometry and structure. 
Distinguishing among the various models is best accomplished by 
making geodetic measurements npt only pear the fault but out to 
distances egual to several times the fault depth. This is where 
the model differences are greatest; these differences will oe 
most readily detected shortly after an earthquake when 
viscoelastic effects are most prono un ced. For a thin channel 
asthenosphere model we .have found that the predicted 
displacements are less than those for a half-space asthenosphere. 
This result is contrary to recently published results based on 
analytical approximations. The deficiencies of the latter work 
result from ignoring material below the asthenosphere and in 
using thickness averaged deformation parameters. Athougn the 
displacements predicted for a thin channel asthenosphere are less 
than those for a half-space asthenosphere, the postseismic strain 
rates at intermediate distances from the fault are greater (in an 
absolute value sense) in the former modeil. 
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CRUSTAL DEFORMATION, THE EARTHQUAKE CICLE, AND MODELS QF 
VISCOELASTIC FLOW IN THE ASTH ENQSPHERE 


INTRODUCTION 

There are Many reasons for the current interest in the 
developient of models of the crustal deformation in active 
seismic zones. Obviously a thorough understanding of the 
deformations that develop prior to an earthquake would contribute 
greatly to predicting a forthcoming event and mitigating the 
hazards associated with it. The analysis of, the deformations 
that accompany the earthguake provide very basic information on 
the rupture process, the size of the event, and the state of 
stre.ss in the lithosphere. Postseismic deformation has been used 
to study relaxation processes within the asthehosphere and the 
coupling between surface slip and aseismic creep at depth.. 

The models that, have been developed to account for the crustal 
deformation of the earthquake cycle fall into two broad 

catagories, the first involving aseismic slip and the second 
viscoelastic flow. The purpose of this paper is to analyze in 
some detail various viscoelastic models of crustal deformation- 
Our objectives are three-fold. First we want to elucidate which 
model parameters are likely to be determinable from present or 
future geodetic measurements. Second we want to compare various 
models to determine whether the different predictions of 
alternative models arise from different physical assumptions or 
different mathematical treatments. Finally we want to determine 
the range of conditions over which substantial viscoelastic 
effects are likely to be present. He are primarily interested in 
the effects cf either vertical or lateral confinement of the low 
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viscosity flew in either the asthenosphere or an .intracrustai 
decoupling lay or separating upper and lower elastic lithosphere 
layers- licdels that focus of these problems have either 
horizontal or lateral variations in viscosity. All of the models 
that we have investigated use linear viscoelastic rheologies and 
all are applied to strike-slip faults. 

The first of the models relevant to the present study was 
developed by Savage and Prescott (1978) for the deformations in 
an elastic lithosphere layer lying over a viscoelastic 
asthenosphere half-space. A conceptually similar, but 
mathematically distinct, model was developed by Spence and 
Turcotte (1979) whose analysis yields similar results. Th.e 
Savage and Prescott model, which we will discuss in more detail 
below, serves as the basic visGoeiastic half-space asthenosphere 
model far our comparisons with other models- These early papers 
pioneered the extern tion of elastic models of the earthquake cycle 
to include viscoelastic flow at depth. Since their publication 
there have been a number of papers in which the vertical and 
horizontal variation of rheological properties has been more 
complicated than an elastic layeo over viscoelastic half-space. 
Prominent among these are the thin asthenosphere model of Leaner 
and Li (1982) and a model with a thin lithosphere near the fault 
and thicker one further away (Yang and Toksoz, 1980) . The latter 
model is called the Varying lithosphere model in this paper. 
These models, along with some variations, will be the subject of 
the comparisons discussed here. Other relevant work has focustd 
on viscoelastic inclusions in subductiop zones (Wahr and Wyss, 
1.972), an intracrustai low- viscosity zone (Turcotte, et.al., 
1983), viscoelasticity both within th,e lower lithosphere and 
asthenosphere (Cohen, 1982) , power law flow within the 

asthenosphere (Helosh and Baefsky, 1983) , and sea anchor effects 
of subducting slabs (Melosh and Fleitout, 1982). All pf the 
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aforementioned models use viscoelasticity as the mechanism for 
time varying deformation rates- Models involving time varying 
aseisais slip are also important, particularly in modeling 
postseisraic rebound. These is, however, a so far unresolved 
ambiguity in distinguishing between elastic slip and viscoelastic 
flow at depth on the basis of surface geodetic measurements. 
Finally it may be important in some locales to take into account 
distributed yielding and the interaction of multiple faults . An 
analysis similar to ours could be made for dip^slip faults by 
building on the modeling work of Bischke (1974), Thatcher and 
Rundle (1979), and Savage (1983).. 


MODELS AND TECHNIQUES 

Tha comparisons that we will present are based on the 

four viscoelastic mo 'els and one elastic ip/tel shown in Figure 1. 
In each case there is a vertical strikes lip fault that extends 
from tha surface down to a depth D. The fault penetrates all or 
part of an elastic layer which we identify with the crust or 
lithosphere- In the simple elastic half-space model this elastic 
layer extends to infinite depth- In the viscoelastic half-space 
model the elastic layer extends from the, surface to a depth H (H 
> Df where a viscoelastic half-space begins.* This low viscosity 
layer is identified with the asthenosphere or an intracrustal 
decoupling layer within the lithosphere. Whereas in Figure lb, 
the low viscosity layer is a half-space, in Figure 1c it has 
finite thickness. The new parameter introduced in this thin 
channel asthenosphere model is the asthenosphere thickness, AH. 
While this model allows for some vertical variation in the 
viscosity structure, other models have allowed for lateral 
variations. One such model, the varying lithosphere model, is 
show in Figure Id. Here the lithosphere is thinner in the 
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vicinity of the fault and thickens farther away* The parameters 
introduced by this model are the lithosphere thickening, A L ai)d 
the half-width of the near- fault, thin lithosphere zone, W. 
Finally Figure 1e is one example of lateral and vertical 
confinement q£ the low viscosity zqne, This inclusion model is 

defined by the half-width panaiieter, w, and the thickness 

parameter, AH* 

The technique we use to compute surface crustal deformations 
com<?s from a finite element adaptation of the procedure used hy 
Savage and Prescott ( 1978) for analysing the repetitive cycle of 
earthquakes on strike-slip faults. The displacements 
attributable to a single episode of strain accumulation agd 
release are decomposed into two contributions. The first 
contribution is due to block motion of the lithosphere at a 
constant velocity, v. The second contribution comes from 
backward fault creep from the surface tq depth, D, at velocity, 
v. The backward creep starts at time zero and continues until 
the earthquake occurrence (recurrence) time,. T. At T an 
earthquake occurs resulting in a relative displacement across the 
fault of AU *= 2vT, strain accumulated prior to the earthquake is 
relieved by the coseismic slip. Considering the combined 
contributions of the two components, the fault is locked at at 
all times before an earthquake as the two terms cancel. At the 
time of the earthquake a sudden relative displacement occurs 
across the fault; subsequently the fault is again locked, A 
sequence of recurring earthquakes is generated by superposing the 
time dependent deformations from individual earthquake cycles 
appropriately shifted in time from one another. The degree bo 
which the pth earthquake cycle in a sequence contributes 
deformation later in the sequence depends on the ratio between 
the recurrence time for earthquakes and the relaxation time of 
the viscoelastic medium. The finite element code that we use in 
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impfeieating the models has been, discussed in previous papers 
(e- g. Cohen, 1982)- We employ a twc dimensional, three degrees 
of freedom, antiplane strain analysis- The rheological 
representation used for the viscoelastic elements is a Maxwell 
fluid in shear and elastic solid in bulk- The Maxwell relaxation 
time, ?•, is defined by r * 2 7) /jlX where 1} is the viscosity 

and yU the rigidity of the viscoelastic material- To facilitate 
comparisons with earlier work, the depth of the coseismxc slip 
plage D will be taken to be the bottpm of the elastic iay«*r, H {D 
= H ia Figures 1b to 1e)- Except in the special cases discussed 
later, the dimensions of the finite element grid ace chosen to be 
large compared to all other dimensions discussed in this paper. 
This assures that the boundary conditions at the sides and the 
bottom of the grid do not effect the results- Explicit 
integration c f the constitutive equations is employed with the 
integration step size chosen to be small compared to the Maxwell 
time of the viscoelastic material- 

It follows from the preaeeding discussion that the 
displacement, u(t), at time t aan be written: 

<*> 

u(t) * vT + u 2 (t) +£ Cu 2 lt + i!5T) - u a (mT)J. (1), 

In this equation the first term is tbq displacement due to the 
block motion and does not contribute to the stress or strain 
field. The second term is the contribution to the crustal 
deformation from the coseismic slip and the aseismic backward 
creep during the current earthquake cycle. The terms under the 
summation sign arise from the viscqelastic response to slip agd 
backward creep during earlier cycles. For an elastic earth each 
of tha difference terms (for a given m) would be exactly zero. 
In addition for an elastic earth u 2 (T) •= 0^ Figure 2 shows how 

the viscoelastic flow in the asthenosphere alters this behavior. 
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Results u co shown tor both naif-space and thin channel 
aathenospheres. Here the individual difference terras are no 
longer zero- The higher order terras in the viscoelastic response 
become of greater relative importance as distance from the fault 
increases. Also tne higher order terras have a somewhat greater 
siynit icanee at most locations in, the thin asthenosphere? model 
than in the half-space case. That the overall amplitude of the 
viscoelastic effect in tne thin channel model is less than in the 
half-space model is also apparent- Summing together all the 
contributions to the displacements ye generate a set of curves 
for the displacements as a function of distance from the fault 
and time between earthquakes as shown la Figure 3. For the 
viscoelastic models, postseismie displacements (e.g. t/T “ 0.2) 
at intermediate distances from the fault not only exceed those in 
the near and far fields but iuso change with a velocity greater 
than the uniform translation rate of the loading blocks- This 
accelerated motion is a common characteristic of the postseismie 
rebound. To raaxe this point more explicit we show in Figure 4 a 
sample plot of velocity, u, versus distance at various times. 
Notice that at t/T =0.1 the velocity at X/D = 2 exceeds that at 
great distances from the fault by over 50 perceut. At some later 
times in the earthquake cycle, however, the velocity is 
considerably slower than that for an elastic earth- 


Let us u ow turn to a more detailed comparison of the thin 
channel anu half-space asthenosphere results. Examples of the 
displacement patterns are compared in Figures 3c and 3d. 
Compaced to the half-space asthenosphere model, the displacements 
in the thin channel model are reduced in amplitude for all times 
t < 1. The dirferences between the two models ace greatest early 
in the cycle, when viscoelastic effects are most pronounced, and 
progress toward zero as the recurrence time is approached- The 
location of the peak in the displacement pattern marks the 
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boundary between a zone of positive shear sfcrAig near the fault 
and ueyative shear strain farther ayay. The width of the zone of 
positive strain increases with time and eventually encompasses 
the entire region shown in the figure* For the thin channel 
model, the location Q f the point of greatest displaceaent is 
displaced faulty ard. This reflects the general confinement of 
the zone of significant deformation (high strain) closer to the 
fault. Toward the end of the cycle, however, the deformation 
field' approaches that for a half-space asthenosphere. Figure 5 
shows an example of the variation ig displacement at X - D as All 
increases. For example with t/T - Or. 2, u/AU * 0.05, 0.10, and 
0.12 for AH/D * 0, 1, and 4 and respectively. The latter value 
is close to the half -space asthenosphere value. 

Since strain rates rather than dispiacemen ts are the 
quantities most readily deduced from most raul tila teration 
measurements we will, henceforth, focus most of our attention on 
a comparison of computed strain rates. Examples of 
representative engineering shear strain rates, t, are shown as a 
function of distance from the fault* X, and at two times in 
Fiyura 6. . The primary observation that can be made about the 

strain rates in the two models is that they are in fairly good 
agreement. This is true despite the fact that the thickness of 
the low viscosity layer in the thin asthenosphere model is fairly 
small, i . e . , AH ~ H. On a closer examination, however, we tigd 
that there are some important second order differences between 
the models. In the near-field, the strain rates after the 
earthquake are somewhat reduced when the asthenosphere thickness 
is decreased. The situation is reversed toward the end of tne 
earthquake cycle. Thus the gear- field strain rates go through a 
greater rang^ of values as the asthenosphere thickness increases. 
Furthermore the minima in the postseismic strain rate patterns 
vary in both magnitude and location with the asthenosphere 
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thickness- No foe ms on these differences in the next two figures 
and discuss those circumstances for which the th'skness of the 
low viscosity zone might be .estimated from comparisons of 
calculations and observations- Figure 7 shows the dependence of 
the near- field (X/D«0.025) postseismic strain rates on the 
astheuosphece thickness AH/D- The normalized strain rate, 
s/(A0/DT) , rises rapidly from its elastic value (when AH -0) of 
0*3 to near the half -space asthenosphere value of just over 1 as 
AH Increases- Even for the very modest value of asthenosphere 
thickuess for which AH » D, e./(AU/DT) -0.9 and for AH/D > 1 the 
strain rate is nearly the half-space value- Thus measurements 
made in the near- field are sensitive to the asthenosphere 
thickuess only over a narrow range of values* If the thickness 
is more than about one fault depth, the shear strain is close to 
tnat obtained with a half-space* Measurements made farther away 
from the fault have a somewhat greater relative sensitive to tne 
asthenosphere thickness. Figure U shows how the magnitude, e ra , 
and location# X m , of the minimum in the postseismic strain rate 
pattern (at a fixed time) vary as a function o£ asthenosphere 
thickness- The magnitude increases (in an absolute value sense) 
as asthenosphere thickness increases until AH/D > 1 then slowly 
decreases toward the half-spaae value of approximately -0,-7. 
Thus foe most values of asthenosphere thickness the minimum is 
deeper than that predicted for the half-space. The position of 
the minimum differs markedly from that for a half-space when the 
asthenosphere thickness is less than about 3D- For example with 
D =20 km., X m = 50 kflu, 65 km*, and 72 km- for Aii = 20 km.,, 40 

km-, id infinity respectively- The position of the minimum for 
the half-space model is almost 50 percent further from the fault 
than it is for AH/D « 1. Note also that the position of the 
minimum moves away from the fault rapidly for AH/D < 1 as elastic 
behavior is approached and the extremum disappears. 
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The results presented above were obtained with the ratio q£ 
the earthquake recurrence time to Maxwell relaxation time equal 
to five, i.e. # T/r * 5. He have rerun a few of our calculations 
with T/r * 20 to see how decreasing the asthenosphere viscosity 

effects these results. As was expected the postsersmic 
displacements are greater when the viscoelastic relaxation is 
more rapid. Therefore the development of the peak an the 
postsaismic displacement pattern at intermediate distances from 
the fault is more pronounced. The postseisaic strain rates 
increase roughly linearly with T/r- Associated nith the increase 
in the postseismic strain rates near the fault is a decrease in 
the rates later on; thus the greatest fraction of the straining 
before an event in an earthquake seguence occurs shortly after 
the proceeding event in the seguence. The strain rate in the 
near-field undergoes an increasing time variation when the 
viscosity of the asthenosphere decreases. Even though the 

postssisiic strain rates are increased with the more rapid 
'"V><J>>9Lastic response, we still figd that the gear- field strain 
rate becomes 90S of the half-space asthenosphere result when AH /D 
> 1-75. This is because the postseismic strain rates increase 
with a reduction in viscosity for bqth the thin channel and half- 
space asthenosphere models. 

One of the most important findings of this study is that the 
effect of having a finite width asthenosphere is much different 
from that predicted by the recent analytic approximations of 
Lehnec and Li (1982), Their model predicts displacements which 
are generally much greater than that for a half-jspaae 

asthenosphere (when X/D < 7) ; whereas, we have found smaller 
displacements when the asthenosphere has a finite width. In 
addition Leaner and Li find that deviations from half-,space 
asthenosphere behavior occur over a much larger range of 
asthenosphere thickness than are indicated by the finite element 
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calculations- For example, with AH/D * 5, they find that the 

displacements at t /£•**<)• 2 ace as much as 40% bigger than that 
predicted by the half-space astheuospere model. on the other 
hand, we find the displacements ace close to that predicted by 
the halt-space ant henosphece model and slighly smaller. After 
performing a number of numerical experiments we now believe that 
we undecstand at least some of the difficulties with the Lehner 
and Li model- Their model is a two-layer approximation with the 
material below the asthenosphere ignored. The mathematical 

*' r 

development q£ the model utilises thickness averaged deformation 
and stress variables and ignores some of the shear stresses 
within the asthenosphere. As the authors point put, the model is 
most suspect near the fault, the region where the greatest 
crustal deformation occurs- The authors also indicate that the 
model may be inaccurate at large distances from the fault- We 
have attempted to mimic the calculatipns of Lehner and Li by 
using their choice of viscosity, rigidity, earthquake recurrence 
time, layer thickness, etc- in a number of calculations. One of 
these calculations is a a two-layer calculation using finite 
elements with the- bottom boundary of the finite element grid held 
fixed for the calculation of u 2 . In a second case, the same 
calculation is performed with the bottom boundary made free- Xu 
a third calculation rfe use out three layer model with elastic 
material extending from the bottom of the asthenosphere to a 
great depth (the choice of the boundary condition at the bottom 
of this grid is immaterial). We have also reproduced the 
analytic calculations and performed calculations using a half- 
space asthenosphere model. (In the latter case several different 
techniques yield substantially the same result. In redoing the 
analytical calculation a sign error was corrected in the 
published equations, but this error has not affected the 

figures.) A comparison of these five calculations is shpwn in 
Figure 10- We see that the analytical model predicts the 
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yr&atest displacements. The smallest displacements are predicted 
by the three layer model and the two layer model with the bottom 
boundary held fixed (these two calculations are in very close 
agreement) , In addition the three layer results are fairly close 
to and a little smaller than the half-space asthenosphere results 
as expected from our earlier analysis. Interestingly the free 
boundary* two-layer model yields displacements which are somewhat 
greater than the half-space asthenosphere displacements. The 
free boundary condition results are the closest results we have 
obtained to the analytical model* but the agreement between the 
two is still not satisfactory. The results obtained with a free 
(fixed) boundary can also be obtained using a three layer model 
with a low (high) rigidity elastic layer below the asthenosphere. 
The validity of the finite element calculations is confirmed not 
only by checking with independent calculations in the half-, space 
asthenosphere limit* AH - cP * but also by verifying the smooth 
transition to elastic behavior as AH approaches zero. on the 
basis of these calculations we conclude that the enhanced 
displacements predicted by the analytical model result from model 
simplifications that are unlikely tq pccur in nature, such as 
ignoring the presence of material below the asthenosphere and 
averaging variables over the fault depth. Despite the fact that 
the displacements are smaller in the thin channel model than in 
the half-space asthenosphere model, we find that the miuima in 
the postseismic strain rate pattern are deeper as indicated 
earlier in Figure 8. 

Before leaving consideration qf the finite width asthenosphere 
model we examine how the halt-width of the strain rate field 
varies with time. lie define the half-width* as the distance 
from the fault at which the st.ta.in rate is ope-half its maximum, 
near-fault value. The width increases with time, an example 
being shown in Figure 11. I»n general the growth in the width of 


ORIGINAL PAGE 63 
OF POOR QUALITY 


►I 


the strain field accelerates with tine following an earthquake as 
the near-field strain rates decrease and the far- field values 
increase- The details of the growth depend op the thickness of 
the asthenosphere. Just after an earthquake the width of the 
strain rate field and its time rate of change are comparable for 
the half-space asthenosphere and the thin layer asthenosphere 
models. As time passes, however, the growth rate for the half- 
space asthenosphere model increases more rapidly. 


The next ornodel we consider, the varyiug lithosphere model, is 
shown in figure Id. For distances from the fault X > W the 
lithosphere thickness is H. At greater distances from the fault 
the thickness is H ♦ AL. As an illustrative case we take W = D 
and AL/D “1.5 (with H/D = 1 and T/r =5). In Figure 12 we once 
again see that the strain rate curves foe this model and the 
half-space asthenosphere model are quite similar but there are 
some important differences. The postseismic near afield strain 
rates are less in the varyiug lithosphere model- Toward the end 
of the earthquake cycle these strain rates are greater- As for 
the spatial pattern, the strain rates at intermediate and far- 
field distances show less variation with distance from the fault 
than do the corresponding strain rates for the half-space model. 
This is because the lithosphere is thicker at these greater 
distances and the deformation pattern takes on some of the 
aspects of <m elastic field. As in the previous case it is 
informative to examine how the magnitude of the near- and 
intermediate-field postseismic strain rates vary with the model 
parameters* Figure 13 shows the variation in the near- field rate 
with width apd thickening. Strain rates within 90% of the half- 
space asthenosphere value are achieved whenever the width of the 
thin lithosphere zone exceeds one-and-a-half times tne fault 
depth, D- Similarly whenever! the distant thickening of tiie 
lithosphere is less than 50%, the strain rates are 90% or greater 
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of the half-space asthenosphere value. As to the location an,d 
magnitude of the postseismic strain rate minima,, the results are 
shown in Figure 14. Four points should be noted. First when the 
width, W/D, is small the magnitude and position of t n depend 
strongly on thickening, AL/D- Second when the width is greater 
than about 31) nearly half-space asthenosphere behavior is 
achieved independent of the value of AL-. Third, when the 
thickening of the lithosphere away from the fault zone exceeds 
three times the fault depth, then 6 m agd X B depend strongly on U 
(for U/D < 3). Fourth, at the small lithosphere thickening value 
of AL = 0.5D, «n/(AU/DT) increases from -.05 to -.71. as W/D 
increases from two tp infinity. The location of the minimum, 
X ra /D, decreases from 4.9 at U/D =0.5 to about 3.5 for U/D egual 
to 2.0 and 3.75 and varies only slightly thereafter. 


The final model we investigate is somewhat of a hybrid of the 
preceeding two models. Here viscoelastic flow is confined to a 
rectangular inclusion located symmetrically about the bottom of 
the rupture tip. The inclusion width, AH, is analogous to the 
layer thickness for the thin layer asthenosphere model; the 
inclusion half-width is w. Figure 15 show,s representative 
calculations of the strain rate for this model with AH/D = 1.5 
and w/D = 1. The viscoelastic features are still present but 
somewhat muted compared to the half-space asthenosphere result. 
Particularly noticeable are the decrease in the postseismic near- 
field strain rates, the poorer development of the postseismic 
intermediate and far-field negative strain rate pattern, and the 
higher strain rates in the near field toward the end of the 
earthquake cycle. These are all effects we have seen in other 
models with finite, dimensions fqr the asthenosphere. Figure 16 
provides information on how the ner:,r-f ield strain rate approaches 
the haLf-sp.ace asthenosphere result as AH and w increase. 
Basically we find that when the thickness and the half-width are 
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gceater than about one or two times the fault depth, tuen the 
near-field strain rates are close to the half-space results. On 
close inspection of the results*, such as those presented in 
Figure 15, ye are again led to the conclusion that the best hope 
for determining the astheno sphere structure cones from making 
very accurate measures at distances up to several times the fault 
depth away flora the fault. 


SUMMARY 

We have examined the strain rates and, to a lesser extent, the 
displacwent patterns associated with several different linear 
viscoelastic models of the earthquake cycle- These models differ 
in the geometric structure assumed for the subsurface 
viscoelastic material. Generally we find that when the 
viscoelastic material has an extent greater thag a few times the 
fault depth, then the crustal deformation pattern looks very much 
like that obtained When the viscoelastic region is a half-space. 

In this case it is difficult to determine information about the 

* 

vertical or lateral structure of the low viscosity zone. On the 
other hand, when the low viscosity zone extends over a smaller 
domain in either a vertical or lateral direction then geodetic 
measurements might be used to , distinguish among models and 
determine the geometric structure of the, asthenosphere. Although 
such a determination will not be unique, the greatest amount of 
information wi.ll come from geodetic surveys made shortly after an 
earthquake (when the viscoelastic effects are greatest) and at 
distances away from the fault up to several times the fault depth 
(where model differences are most important). Among the features 
predicted by the finite element viscoelastic calculations is that 
the confinenent of viscoelastic flow to a finite thickness 
channel reduces the displacements over that expected for a 
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viscoelastic half-space. Indications to the contrary .based on 
analytic work; are the result of ignoring the presence of high 
viscosity Material beloy the viscoelastic channel and usiqg 
thickness averaged variables Mhich are inapplicable near the 
fault. 
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Figure Is Models used in earthquake cycle simulations: 

a. Elastic half-space model 

b. Viscoelastic nal£~spac$ asthenosphere model 

c. Viscoelastic thin-channel asthenosphere model 

d. Viscoelastic varying lithosphere model 
e« Viscoelastic inclusion model. 


Figura 2: Superposition of multicycle displacement fields for 

thin channel and half space asthenosphere models. 


Figure 3: Displacements versus distance from the fault at various 
times for several different models: 
a. Elastic model 

Jb. • Viscoelastic half-space asthenosphere model 
c. Viscoelastic thin channel asthenosphere model. 


Figura 4: Velocity versus distance from the fault at various 

times for half-space asthenosphere model- 


Figure 5: Displacement at x/D = 1 and at various times versus 

thickness of asthenosphere channel. 


Figure 6: Strain rate versus distance from the fault at various 

times for half-space and thin channel asthenosphere models. 


Figure 7: Near -field, postseismic strain rate versus 

asthenosphere thickness for thin channel asthenosphere model. 


Figure 8; Magnitude and location of postseismic strain rate 
minimum for thin channel asthonosph^re model. 
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Figure 9; Comparison of displacements and strain rates for 
viscelastic models with different relaxation times. 

Figure 10: Half-width of strain rate field versus time for half- 

space and thin channel asthenosphere models. 

Flyura 11: Comparison of patterns of displacement versus distance 
from the fault for thin channel analytic model and finite element 
models employing various boundary conditions (b. c.) for the 
calculation of u 2 . 

Figure 12: Strain rate versus distance from the fault at various 

times for half-space asthenosphere and varying lithosphere 
models. 

Figure 13: Hear~.fi eld postseismic strain rate dependence on the 

width of the thin lithosphere zone and magnitude of the 
lithosphere thickening. 

Figure 14: Magnitude and position of postseismic strain rate 

minimum for varying lithosphere model. 

Figure 15: Strain rate versus distance from the fault at various 

times for the half-space asthenosphere and inclusion models. 

Figure 16: Near-field postseismic strain, rate dependence on 

inclusion width and thickness. 
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Figure 1: Mpdels used in earthquake cycle simulations: 

a. Elastic half-space model 

b. Viscoelastic half-space asthenosphere model 

c. Viscoelastic tbiq-chan^el asthenosphere model 

d. Viscoelastic varyii^g lithosphere model 

e. Viscoelastic inclusion model. 
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Figure 2z Superposition of multicycle displacement fields for 
thin channel and half space asthenosphere models. 
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Figure 6: Strain rate versus distance fro* the fault at various 

tijn^s for half-space and thin channel astheqosphere Models. 
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Figure 9: Comparison of displacements and strain rates for 

viscelastic models with different relaxation times. 
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Figure 10s Half-vidth of strain rate field versus tine for half- 
space and thia channel astheaosphere models. 


30 



ORIGINAL PAGE 1ST 
OF POOR QUALITY 


U(m) 



Figure 11s comparison of patterns of displacement versus distance 
from the fault for thin channel analytic model and finite element 
models employing various houndary conditions (b.c.) for the 
calculation of u 2 . 
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Figure *2: Strain rate versus distance iron the fault at various 
times for half-space as the nosphere and varying lithosphere 

models q 
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Figure 13; Near-field postseisaic strain rate dependence on tli,e 
width of the thin lithospiiere zone and Magnitude of 
lithosphere thickening. 
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Figure 34i Magnitude and position o£ postseisaic strain rato 
aimiaua for varying lithosphere model- 
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Figure 16: Near-field post seismic strain, rate dependence on 

inclusion width and thickness. 
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